Mach cone induced by 7-triggered jets in high-energy heavy-ion colhsions 
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Medium excitation by jet shower propagation inside a quark-gluon plasma is studied witiiin a 
linear Boltzmann transport and a multiphase transport model. Contrary to the naive expectation, 
it is the deflection of both the jet shower and the Mach-cone-like excitation in an expanding medium 
that is found to gives rise to a double-peak azimuthal particle distribution with respect to the initial 
jet direction. Such deflection is the strongest for hadron-triggered jets which are often produced 
close to the surface of dense medium due to trigger-bias and travel against or tangential to the radial 
flow. Without such trigger bias, the effect of deflection on 7-jet showers and their medium excitation 
is weaker. Comparative study of hadron and 7-triggered particle correlations can therefore reveal 
the dynamics of jet-induced medium excitation in high-energy heavy-ion collisions. 
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Strong jet quenching has been observed in experiments 
[11-3] at the Relativistic Heavy-ion Collider (RHIC) as a 
consequence of jet quenching or parton energy loss in 
high-energy heavy- ion collisions (H . The energy and mo- 
mentum lost by a propagating parton will be carried by 
radiated gluons and recoiled medium partons vifhich in 
turn will go through further interaction and eventually 
lead to collective medium excitation such as supersonic 
waves or Mach cones [1,01 ■ Indeed, Mach cones have been 
found in the solutions of both hydrodynamic response 
8|-ll2| and linearized Einstein equations in string theory 
I3I Il4| excited by a propagating jet. Such collective 
excitation by a propagating jet is expected to be respon- 
sible for the observed conic back-to-back (b2b) azimuthal 
dihadron [isj, [l^ and trihadron correlations [l3| with a 
maximum opening angle of A(j) w 1 (rad) relative to the 
backside of a triggered high-p^ hadron. However, hadron 
spectra from the freeze-out of the Mach cone in both hy- 
drodynamics with realistic energy-momentum deposition 
by jets [1, [1, [l0| and string calculations in the hydrody- 
namic regime |14{ fail to reproduce the observed conic 
azimuthal correlations. Such correlations on the other 
hand are observed in a multiphase transport (AMPT) 
Monte Carlo simulations [l^ which could come from jet- 
induced wakes that are deflected by a radially expanding 
medium [lOj . 

Dihadrons with a high-p^ trigger are mostly domi- 
nated by b2b jets that are produced close to the sur- 
face of the dense matter [l^ with the awayside jets often 
traveling against or tangential to the radial flow. De- 
flection of these jet showers and associated Mach cones 
by the radial flow can lead to double-peaked hadron az- 
imuthal correlations. On the other hand, high-py 7's 
are produced throughout the volume of the dense mat- 
ter [201 . The effect of deflection should be reduced for 
7-triggered jet showers after averaging over all possible 
production positions and propagation direction, leading 
to a weaker double-hump 7-hadron correlation as com- 
pared to dihadron correlation. 

In this Letter, we will study medium excitation by a 



propagating jet shower usin g b oth a linear Boltzmann 
transport and AMPT model [2l|. We will illustrate that 
while a Mach-cone-like excitation by a propagating jet in 
a uniform medium cannot give rise to a conic distribu- 
tion of the final partons, deflection of the jet shower and 
the Mach-cone-like excitation in an expanding medium 
will result in a double-peaked azimuthal distribution as 
observed in dihadron measurements [H, [ll|. Because 
of the different geometric distributions and propagation 
direction of the initial produced jets, we will illustrate 
that 7-hadron and dihadron azimuthal correlation will 
be quantitatively different, depending on the value of 
jet-medium cross section. We therefore propose to use 
comparative study of 7-hadron and dihadron azimuthal 
correlations to shed light on the dynamics of jet-induced 
Mach-cone-like excitation in high-energy heavy-ion colli- 
sions. It can also disentangle other mechanisms such as 
triangular flow and hot spots [1^ that contribute to 
dihadron but not 7-hadron azimuthal correlation, though 
they in general enhance jet-medium interaction and the 
resulting medium modification of 7-hadron and dihadron 
azimuthal correlations. 

We first study the jet shower propagation and medium 
excitation through a linearized Monte Carlo simulation 
of the Boltzmann transport equation 

Pi ■ dflipi) = - J dp2dp3dp4iflf2 - hfi) |A/l2^34|^ 

X [2TT)H\pi+p2-p:i-Pi), (1) 

including only elastic 1 + 2 — > 3 + 4 processes as 
given by the matrix elements Mi2_>.34, where dpi = 
d^p^/[2E,{2T:)% fi = l/(eP'-"/^ ± 1) (« = 2,4) are ther- 
mal parton phase-space densities in a medium with lo- 
cal temperature T and flow velocity u = (1, v)/\/l — , 
fi = {27r)^ (p ~ Pi) 6^ (x — Xi — Vit) [i ~ 1,3) are the 
jet shower parton phase-space densities before and af- 
ter scattering, and we neglect the quantum statistics in 
the final state of the scattering. We will consider quark 
propagation in a thermal medium and assume small an- 
gle approximation of the elastic scattering amplitude 
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\Mi2^3a\ = Cg^{s'^ +u'^) / {t+ fi"^)"^ with a screened gluon 
propagator, where s, t and u are Mandelstam variables, 
C=l (9/4) is the color factor for quark-gluon (gluon- 
gluon) scattering and fi is the screening mass which 
we consider here as a constant cut-off of small angle 
scattering. The corresponding elastic cross section is 
da/dt = |Afi2-^34|^ /Wtts^. 

We assume partons follow classical trajectories be- 
tween scatterings and keep track of the propagation of 
both jet shower partons and the recoiled medium parton 
after each scattering. Both the jet and recoiled medium 
parton are allowed to scatter again with other medium 
partons as they further propagate. We determine the 
probability of jet-medium scattering in each time step 
Atj according to an exponential distribution 
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P, = 1 - cxp[- Axj-uaip{xj,tj)], 



(2) 



where cr^ is the scattering cross section between jet and 
medium partons, the sum is over the time steps along 
the classical trajectory since the last scattering and p is 
the local medium parton density. For each scattering, we 
keep record of both the leading jet parton (pa) and the 
recoil medium parton (^4). Thermal partons with the 
initial momentum (^2) will be subtracted from final par- 
ton phase-space density to account for the back-reaction 
in the Boltzmann transport equation. The net parton 
phase-space density Sf{p) averaged over many events will 
be the medium excitation by the propagating jet. We call 
such Monte Carlo simulation a linearized Boltzmann jet 
transport since we neglect scatterings between recoiled 
medium partons. This is a good approximation to the 
full Boltzmann jet transport as long as the medium ex- 
citation remains relatively small 5f{p) <C /(p). 

Within such a linearized Boltzmann jet transport one 
can study not only parton energy loss but also the evolu- 
tion of the medium excitation induced by the propagating 
jet. Shown in the upper panel of Fig. [l]arc the contour 
plots of the energy density rdE/drdz of the medium ex- 
citation induced by a quark with initial energy E ~ 20 
GeV at different times after the start of the propagation 
along z direction in a uniform medium with temperature 
T ^ 300 MeV. We used = 0.5 and fi = 1.0 GeV. One 
can clearly see the formation of Mach-cone-like medium 
excitation at later times and a diffusion wake behind the 
jet which has negative energy density due to depletion 
of partons in this region by the propagating jet, simi- 
lar to results of the linearized hydrodynamic response 
S S, [m, d and string theory [Tj, EM calculations. 

We can also calculate parton spectra at different times, 
corresponding to equal-time freeze-out in hydrodynamic 
studies. Shown in the lower panel of Fig. [T]are the az- 
imuthal distributions of partons from the medium exci- 
tation at different times. At early times, the azimuthal 
distribution is dominated by recoil thermal partons from 
the primary jet-medium scattering, which is given by 
the collision kernel in the Boltzmann equation. It has a 
double- hump feature in the azimuthal distribution for low 




FIG. 1: (Color online) Contour plot of energy density 
rdE/drdz of medium excitation induced by a quark with ini- 
tial energy E = 20 GeV that propagates in the z direction in a 
uniform gluonic medium (upper panel) and the corresponding 
azimuthal parton distributions (lower panel). 



energy partons with opening angles determined by the 
value of temperature and screening mass in the scatter- 
ing matrix element. As a parton propagates through the 
medium, such recoil thermal partons alway accompany it 
at the neck of the Mach-cone-like excitation. This is very 
similar to the nonequilibrium part of jet-induced medium 
excitation in the string theory study 124 1 . 

At later times, the recoil medium partons from the 
primary jet-medium interaction will scatter with other 
thermal partons as they further propagate through the 
medium, causing diffusion of the wake front. As a conse- 
quence, low energy partons from the body of the Mach- 
cone-like excitation overwhelm that from the neck and 
the final azimuthal distribution has only a broad single 
peak along the direction of the propagating jet. 

The hot matter created in high-energy heavy-ion col- 
lisions has a finite initial transverse size as given by that 
of two colliding nuclei. Because of the tremendous initial 
pressure, the hot matter will experience rapid transverse 
expansion and develop strong radial flow. The nonunifor- 
mity over a finite transverse size and strong radial flow 
of the hot medium should influence the propagation of 
jet showers and their medium excitation. 

To take into account the dynamical evolution of the hot 
matter in our linearized Boltzmann jet transport model, 
we will use the numerical results from a (3-t-l)D ideal hy- 
drodynamical calculation (25| for local temperature and 
flow velocities, from which we generate thermal momen- 
tum p2 and Pi for each jet shower and medium parton 
interaction in Eq. ([1]). We also use HIJING Monte Carlo 
model [2^ to provide the initial jet shower parton dis- 
tribution which consists of multiple partons from the 
final-state radiation in 7-jet events. Shown in the upper 
panel of Fig. [5] are the contour plots of the energy den- 
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sity in both transverse (x-y) and the beam (a;-^) plane 
excited by a 7-triggered jet shower with initial position 
{x, y, z) ~ (—4, 0, 0) fm and energy Ej. = 20 GeV in cen- 
tral Au-\-Au collisions at ^/s = 200 GeV/n. As compared 
to the case of a uniform medium, the shape of the medium 
excitation is distorted considerably by the transverse and 
longitudinal flow of the expanding medium. The distor- 
tion depends on the direction of the jet propagation rel- 
ative to the flow. 

y[fm] z[fm] 




FIG. 2: (Color online) (upper panel) Contour plot in the 
transverse {x-y) and beam (x-z) plane of energy density ex- 
cited by a quark jet shower with Et = 20 GeV and initial 
position at {x,y,z) = (—4,0,0) fm that travels toward the 
center of the expanding medium as given by ideal hydrody- 
namics [2^ for central Au + Au collisions at the RHIC energy. 
The azimuthal distribution of medium and jet shower partons 
when the jet shower travels against (lower left) and perpen- 
dicular (lower right) to the transverse flow. 

The azimuthal distribution of partons from both the 
jet shower and jet-induced medium excitation is also dis- 
torted by the transverse flow and the nonuniformity of 
the dense medium. To illustrate the influence of trans- 
verse flow and the density gradient, we show in the lower 
panel of Fig. [2] the azimuthal distribution of jet shower 
and medium partons from a 7-triggered jet that is pro- 
duced at an initial position (x, y, z) = (—4, 0, 0) fm away 
from the center of the dense medium and propagating 
against (lower left) and perpendicular (lower right) to 
the radial flow, respectively. For a tangentially propa- 
gating jet shower (lower right), low px partons from the 
jet shower and Mach-cone-likc excitation arc clearly de- 
flected by both the density gradience (we verify this by 
setting the transverse flow velocity to zero) and the ra- 
dial flow, giving rise to the azimuthal distributions that 
peak at an angle away from the initial jet direction. For 
jet showers that travel against the radial flow (lower left), 
the same deflection essentially splits the azimuthal dis- 
tribution of low pt partons to become a double-peaked 
one. Such deflection in an expanding system is also ob- 



served in an ideal hydrodynamical study of Mach cone 
propagation and it will give rise to both the diag- 
onal (tangential jet showers) and off-diagonal (split jet 
shower) part of the 3 particle correlation if analyzed as 
in the experimental study [l7| . The magnitude of the 
conic correlations depends on the parton cross section 
while the opening (deflection) angle also depends on the 
radial flow velocity. 

To obtain the final jet-induced medium parton distri- 
bution in 7-jet events in heavy-ion collisions, one should 
average over the propagation direction and the initial 
production positions that are distributed in the trans- 
verse plane according to the binary collisions with the 
nuclear geometry. The averaged azimuthal distribution 
is found to have a broadened single peak in the jet direc- 
tion for small values of the parton cross section. The de- 
flection is expected to give a double-peak azimuthal par- 
ton distribution for sufficiently large values of the parton 
cross section. However, within the linearized Boltzmann 
jet transport model, larger parton cross sections will lead 
to larger amplitudes of the medium excitation and even- 
tually the linear Boltzmann transport breaks down which 
neglects interaction among jet-excited medium partons. 
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FIG. 3: (Color online) Dihadron (open symbol) and 7-hadron 
(filled symbol) azimuthal correlation from AMPT [2l| model 
calculation with different values of parton cross section. 

To go beyond the linear Boltzmann transport, we use 
the AMPT Monte Carlo model to simulate hadron and 
7-triggered jet events and study jet-induced medium ex- 
citation. AMPT [2l| is essentially a full Boltzmann par- 
ton and hadron transport model with only elastic par- 
ton collisions and initial conditions given by the HIJING 
model ^26]. Shown in Fig. |3] are azimuthal distributions 
of hadrons (with p!^^'^° = 1 — 2 GeV/c) associated with 
both a high-pT trigger 7 (closed symbols with solid lines) 
and hadron (open symbols with dashed lines) in cen- 
tral Au + Au collisions at the RHIC energy with dif- 
ferent values of the parton cross section. The hadron 
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FIG. 4: (Color online) 7-hadron (filled symbol) and 7-parton 
(open symbol) azimuthal correlation from AMPT [23| model 
calculation with different values of parton cross section. Re- 
sults are scaled with different factors for clear presentation 

correlations develop a double-peak feature in the oppo- 
site direction of the trigger due to the deflection of jet 
shower and Mach-cone-like excitation by the transverse 
flow and density gradient as one increases the value of 
parton cross section. The opening angle between two 
peaks increases with the value of the parton cross sec- 
tion. At very large values of the parton cross section, one 
should approach the hydrodynamic limit with an open- 
ing A(f> sa 1 (rad) for the double-peak structure. Most 
importantly, one can observe that the amplitudes of the 
double-peaks in dihadron correlations are much bigger 
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than that of 7-hadron correlations for a given value of 
parton cross section, verifying our argument that the de- 
flection of jet shower and Mach-cone-like excitation by 
the radial flow has a stronger effect on the azimuthal di- 
hadron correlation than 7-hadron correlation because of 
the difference in the trigger bias on the initial jet pro- 
duction position and propagation direction. Therefore 
comparative study of dihadron and 7-hadron azimuthal 
correlations can shed light on the dynamics of jet propa- 
gation, medium excitation, the strength of medium par- 
ton interaction and other effects such as hot spots [l^l 
and triangle flow (22j in high-energy heavy-ion collisions. 

A hadronization mechanism such as parton coalescence 
employed in AMPT model can quantitatively influence 
the final hadron correlation [l^j as compared to the par- 
ton correlation before hadronization. As shown in Fig. |4] 
such an effect, as manifested in the differences between 
7-hadron and 7-parton correlation before hadronization, 
is rather small in AMPT. We have not included inelastic 
parton interaction in this study but expect the results 
remain qualitatively the same. However, for quantitative 
theoretical predictions and comparison to experimental 
data, one should include inelastic parton interaction and 
hadronic interaction in future studies. In addition, the 
finite formation time for the jet parton shower, which 
is neglected in this study, will delay the jet shower and 
medium parton interaction and quantitatively influence 
the final jet-induced 7-hadron and dihadron correlations. 
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